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We propose a simple way to parameterize the gap function in iron pnictides. The key idea 
is to use orbital representation, not band representation, and to assume real-space short-range 
pairing. Our parameterization reproduces fairly well the structure of gap function obtained in 
microscopic calculation. At the same time the present parameterization is simple enough to 
obtain an intuitive picture and to develop a phenomenological theory. We also discuss simpli- 
fication of the treatment of the superconducting state. 
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What is new in the iron pnictide 1 ) as a high-T c super- 
conducting material is its multi-orbital nature. Several 
microscopic calculations have been carried out to show 
that the five 3d orbitals of Fe atoms are entangled 2 ' 3 ) 
and that the superconducting gap function on the multi- 
ple Fermi surfaces are very complicated. 4 9 ) Therefore, it 
is desirable to construct a description of the gap function, 
which is simple enough to obtain an intuitive picture for 
developing a phenomenological theory, but also which is 
powerful enough to capture its complicated multi-orbital 
nature. In fact, the gap functions of several kinds of iron- 
pnictide superconductors are not universal, i.e, nodeless 
in some materials and nodal in other materials. 10 18 ) Fur- 
thermore, there are some controversial results on the ex- 
istence of inter-band sign reversal of the gap function. 
For instance, although quasi-particle interference mea- 
surements imply the existence of the inter-band sign re- 
versal, 19 ) robustness of T c against impurity scattering 
suggests that there is no sign change between bands. 20 ' 21 ) 
Having these experimental situations in mind, we aim to 
construct a simple description of the gap function that 
enables us to discuss such complexity in a simple man- 
ner. We use orbital representation instead of band rep- 
resentation, and assume real-space short-range pairing. 
Our parameterization reproduces very well the structure 
of gap function obtained in microscopic RPA calcula- 
tion. We also discuss simplification of the gap function 
for studying the superconducting state of iron pnictides. 

We analyze the structure of the gap function based on 
the multi-orbital Hubbard model proposed by Kuroki et 
al. 3 ^ that is downfolded from first-principle calculation. 
In this downfolding scheme, five bands around the Fermi 
energy are kept. The obtained five basis wave functions 
have the symmetry of Fe-3d orbitals, i.e., d 3z 2_ r 2, d zxi 
d yzi d x 2_ y 2 and d xy . One thing to be noted here is that 
the extended Brillouin zone is used, i.e, only one Fe atom 
is contained in a unit cell. We use the Hamiltonian H = 
Ho + Hi with 

^0 = t ^Mi,3 C \aa C 3bcj = S ab{k)c ] ka(J C kh(J , (1) 

i,j,cr,ab ka,ab 



where t a ,6;i,j is that obtained in the downfolding and 
Hi is the standard onsite multi-orbital interaction with 
(U,U' J#, J). 22 ' 23 ) In eq. (1), indices a and b run 
through to 4, which correspond to d 3z 2_ r 2, d zxi d yzi 
d x 2_ y 2 and d xy orbitals respectively. Details of the hop- 
ping integrals are different from system to system. In this 
paper, we mainly use the model with hopping integrals 
shown in Table I of Kuroki et aZ., 24 ) in which the hopping 
integrals up to the fifth-nearest neighbors are kept and 
the three dimensionality is neglected. We call this model 
as typical-1111 model in the following. 25 ) Full models for 
LaFeAsO, LaFePO and NdFeAsO are also used. 

Before analyzing the gap functions, we discuss the 
orbital- and band-representations in the multi-orbital 
Hubbard model. The Hamiltonian, eq. (1), is written in 
the orbital representation and Hq can be diagonalized 
by a momentum-dependent unitary transformation. We 
call the transformed representation as band representa- 
tion, in which the free Green's function, Gn(fc) (5x5 ma- 
trix form, k = (i(j n ,k)), is diagonalized. However, since 
the momentum dependence of the unitary transforma- 
tion is very severe in iron pnictides, the orbital character 
strongly depends on the Fermi surface positions. This 
causes the complexity when we use the band representa- 
tion. In contrast, we find that the orbital representation 
is often suitable to discuss physics in real-space view. 
Although the Green's function, Go(fc), in this represen- 
tation has off-diagonal components connecting different 
orbitals, it does not make serious problems when we use 
the matrix form of Dyson-Gor'kov equation, 

((G^k))- 1 A k \ (6(k) F(k) \ - 
V At -(t&>(-k))- 1 ) [P( k )-tG(-k)) ' 

which is in the present case, 10x10 matrix. Note that we 
only consider the singlet channel and neglect the normal 
self-energy and frequency dependence of the gap function 
in eq. (2). 

Now we begin to investigate the structure of A^. The 
key idea is to use the orbital representation and to con- 
sider in both real and momentum space. From now on, 
Afc is always in the orbital representation. First, we con- 
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sider the diagonal elements of A&, i.e., the intra-orbital 
pairing in momentum space. Limiting our consideration 
to the s-wave channel, what is important is the gap value 
at (0,0),(7r,0) and (7r,7r) in the extended Brillouin zone, 
since the Fermi surfaces of iron pnictides are small and 
enclose those points. We denote the gap values on the 
orbital i at (0,0), (tt,0) and (tt,tt) as A r , Af and Af, 
respectively. For later use, we also define the difference of 
the gap value at (±7r/2, tt/2) as A+ — A~ (see Fig. 1(a)). 



ble pairing between the different orbitals (Fig. 2), we 
obtain 




(b) c-S, 
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Fig. 1. (Color online) (a) Schematic picture of the definition of 
Af, Af, AV' and Af. (b) Real space picture of the Aa(k). 



Now we move on to the real space picture. In order 
to make a simple description of the gap functions, we 
assume short-range pairings in the orbital representation. 
Then, Aa(k) becomes 

Au(k) = cii + 2bi (cos k x + cos k y ) + 4q cos k x cos k y 

+4^ sin k x sin k y , (3) 

where we have introduced three parameters, ai, bi and 
Ci representing on-site (a^), nearest-neighbor (bi), and 
next-nearest-neighbor (q) pairings. Anisotropy parame- 
ter Si is introduced for the d zx and d yz orbitals (i =1, 
2), since they do not have 4- fold symmetry, i.e., we have 
Si = — 62 = S and Sq = S3 = £4 = 0. Meanings of these 
parameters are described in Fig. 1(b). Relation between 
the parameters (Af , Af 1 , Af , Af) and (ai, bi,ci,Si) can 
be simply written as 

ai = (Af + 2Af + AD/4, bi = (Af - AD/8, 

a = (Af - 2Af + AD/16, ^ = (A" - A+)/8. (4) 

Note that the symmetry of 3d orbitals give a\ = a 2 , 
bi = b 2 , ci=c 2 . 

Next, we consider the off-diagonal elements of A&, i.e., 
inter-orbital pairings. Assuming only the shortest possi- 
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Fig. 2. (Color online) Up left: Symmetries of the basis wave func- 
tions. Note that the coordinate of the original unit cell is used 
to describe the basis wave functions and they are rotated by 
45° from the usual definitions of the d-orbits on the square lat- 
tice. Others: Schematic real space pictures for the off-diagonal 
elements of A(fe). Red circle represents the As sites. 



Aoi(fc) 


= 2ip i(sinfc a; - 


sin k y ), 


(5a) 


A 02 (fe) 


= -2ip 01 (smk x 


+ sin fcjj, 


(5b) 


A 03 (k) 


= 2d 03 (cos(k x + 


fc y ) - cos(^ - k y )), 


(5c) 


A 04 (fc) 


= 2^04 (cos k x — 


cos & y ), 


(5d) 


Ai 2 (fe) 


— 2di2(cos k x — 


cos fc^), 


(5e) 


Ai 3 (fe) 


= 2ip 13 (smk x - 


sin k y ), 


(5f) 


A 23 (fe) 


= 2ip 13 (smk x + 


sin fey), 


(5g) 


Au(fe) 


= 2ip 14 (sink x + 


sin k y ), 


(5h) 


A 24 (fe) 


= —2ipi4(sm k x 


— sin k y ). 


(5i) 



Here the symmetries of orbitals are important for the sign 
of pairings as shown in Fig. 2. Furthermore, the matrix 
elements concerning d zx and d yz orbitals (i =1, 2) are 
determined by taking account of As atoms which exist 
above and below the Fe plane, because the d zx /d yz or- 
bitals are odd in z-direction and they do not have matrix 
elements with others without As. Characteristic feature 
found in eqs. (5) is that the off-diagonal elements have 
p-wave or d-wave symmetry. This is a direct consequence 
of the symmetry of the basis wave functions and the fact 
that the diagonal elements have s-wave symmetry. In the 
case where the diagonal elements have d-wave symmetry, 
we can do the similar analysis. 

In order to determine the parameters in eqs. (3) and 
(5), we calculate the gap function A^ PA , by solving lin- 
earized Eliashberg equation with the effective interac- 
tion obtained within RPA. In this calculation, we fol- 
low the formalisms in Kuroki et al. 3 ^ We use the filling 
n = 6.1, temperature T = 0.02 eV and Coulomb parame- 
ters (U,U',J H ,J) = (1.0,0.6,0.2,0.2) eV. The Brillouin 
zone is divided into 32x32 meshes and 512 Matsubara 
frequencies are used. Then, using the gap function A^ PA , 
we obtain (A r , Af 1 , A r , Af) from which the parameters 
(ai, bi,Ci,Si) are determined through eqs. (4). Parameters 
Pij and dij are determined from Af J PA (k) (i ^ j) simi- 
larly. The obtained results for a typical-1111 model are 
shown in Table. I. 

In Fig. 3, we compare the present simple description 
An(fe), A 33 (k) and Ai 3 (fc) with those in A£ PA deter- 
mined microscopically. We can see that our fitting is in 
good agreement with A^ PA . Although some of the off- 
diagonal elements in A(fc) do not show a good coinci- 
dence, their magnitudes are small so that it does not lead 
to a serious problem. It is rather surprising that the gap 
function can be expressed only with the small number 
of parameters representing the real-space pairing up to 
next-nearest-neighbor. In contrast, when the gap func- 
tions are re-expressed in the band representation, they 
become very complicated for getting an intuitive picture. 
In Fig. 3(d), we also plot the density of states (DOS) 
calculated with the fitted gap functions. The result shows 
nodeless two-gap behavior where the size of the smaller 
gap is about half of the larger gap. 

We have checked the present parameterization for var- 
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Table I. Parameters obtained in RPA calculation. Typical 6.1, typical 6.3 are the results for typical-1111 model with the filling n =6.1 
and 6.3 respectively. The result for LaFePO is obtained with the filling n =6.1. 





ai 




ci 


S 


CL3 


b 3 


C3 


di2 


P13 


typical 6.1 


0.000746 


-0.00382 


-0.00155 


0.000751 


0.0213 


0.0111 


-0.00474 


-0.000536 


0.00446 


typical 6.3 


0.00185 


-0.000598 


-0.00517 


0.000895 


-0.00541 


0.0117 


0.000315 


-0.000248 


0.00642 


LaFePO 


0.00596 


-0.00301 


-0.00107 


0.000504 


-0.00716 


0.00235 


-0.000441 


-0.00220 


0.00166 




Fig. 3. (Color online) (a-c) Calculated An(fc), A33(fc) and 
Ai3(fe) in RPA and fitted ones for typical-1111 model, (d) 
DOS for typical-1111 model calculated with 2048x2048 k-space 
meshes. 



ious cases to find that the fitting is very good for most of 
the cases. Exception is the case when we have too strong 
spin fluctuation on d x 2_ y 2 orbitals, for example in the 
model for NdFeAsO. In this case, Ass(k) has strong mo- 
mentum dependence and eq. (3) is not enough to repro- 
duce its strong momentum dependence. However, it is 
well known that the RPA estimation of the effective in- 
teraction gives stronger momentum dependence of A(fc) 
compared with other approximations. Therefore, we ex- 
pect that the actual momentum dependence of Ass(k) in 
this case is much milder and the present parameteriza- 
tion works as well. In addition, the damping effect caused 
by impurities and the strong-correlation effects neglected 
in RPA generally give smoother gap functions. 

In the following, we consider the possible simplifica- 
tion of the gap functions. For the typical-1111 model 
and the models for LaFeAsO and NdFeAsO, 25 ^ we ex- 
pect that only d zx /d yz and d x i_ y i orbitals are necessary 
since they contribute most of the DOS near the Fermi 
energy. In this case, we need only parameters (oi, &i, 
ci, 5, 03, 63, C3, di2, P13). We find that DOS calculated 
with keeping these parameters and setting others to zero 
can reproduce the DOS in Fig. 3(d). Furthermore, even 
if we assume d\2 = P12 = 0, the obtained DOS repro- 
duces most of the features in Fig. 3(d), although there 
is a slightly different structure for the smaller gap. Since 
the resultant gap function is so simple, it is very useful to 
obtain an intuitive picture. Note that the above simplifi- 
cation can not be applied when d 3z 2_ r 2 or d xy contribute 
to the DOS around Fermi energy, for example, in the case 
of LaFePO. 

Next, we study the the difference between the nodal 
and nodeless s-wave gap in terms of the present param- 
eterization. The existence of the nodal s-wave is theo- 
retically predicted for the model for LaFePO. 3 ) Thus, 



we carried out similar calculation using the model of 
LaFePO. Here we use three-dimensional model and di- 
vide the Brillouin zone into 32x32x4 meshes. We find 
that the most significant difference from the typical-1111 
model appears in A33(fe), which is shown in Fig 4 (see 
also Table. I). Comparing Figs. 4 and 3, we can see that 
A33(fe) for LaFePO model has weaker momentum depen- 
dence than that for the typical-1111 model. In terms of 
our parameterization, cs is smaller in LaFePO. Due to 
this difference, A33(fc) at (71", 0) is positive in the typical- 
1111 model and negative in LaFePO. Note that An(fe) 
at (7r, 0) is positive for both models. Figure 4(b) shows 
the obtained DOS for LaFePO. We can see that (1) the 
partial DOS for the d zx /d yz and d x i_ y i orbitals show 
nodal behavior and (2) there is large remaining DOS in 
the gap around 20 meV since only a very small gap (< 
2 meV) appears in the d 3z 2_ r 2 orbital. 




Fig. 4. (Color online) (a) Calculated A33(fc) and fitted one for 
LaFePO model, (b) DOS for LaFePO model calculated with 
512x512x32 k-space meshes. 



To obtain a clear view on how the node evolves, we cal- 
culate DOS with a hypothetical gap function, with only 
(ai, 03) = (A, — A). This means that we consider the to- 
tally momentum-independent s-wave gap in the orbital 
representation but with different sign on d zx /d yz orbitals 
and d x 2_ y 2 orbitals. Interestingly, obtained DOS shows 
clear nodal behavior (Fig. 5(a)) although we input the 
gap function without any momentum dependence. This 
is because the orbital character varies along the posi- 
tions of the Fermi surface. Some part of the Fermi sur- 
face enclosing (71- , 0) has d zx /d yz character and the gap 
function is +A, while other part of that Fermi surface 
has d x 2_ y 2 character and the gap function is — A. In this 
case, nodes appear on the Fermi surface. Although the 
gap with (01,03) = (A, —A) is not the result of micro- 
scopic calculation, this argument captures the nature of 
the node obtained for LaFePO model. For comparison, 
we show the result with (01,03) = (A, A) in Fig. 5(b). 
It shows clear full-gap feature and no nodal behavior as 
expected. 

Finally we investigate the doping dependence of the 
gap function within RPA calculation concentrating on 
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Fig. 5. (Color online) DOS for the hypothetical gap function 
of (ai,a 3 ) = (A, A) (a) and (A, -A) (b) calculated with 
2048x2048 k-space meshes and A =40 meV. 



the electron-doping side. We use the typical-1111 model 
and treat the doping as the rigid band shift. As we can 
see from Fig. 6, Ass(k) becomes flatter and flatter with 
doping. Figure 6(d) shows <23, bs and cs in Ass(k) as 
a function of doping. The most important behavior is 
that of C3. With doping the absolute value of C3 be- 
comes smaller and C3 disappears around n ~ 6.25. Ac- 
tually, the strongest superconducting instability channel 
changes from s-wave to d-wave around n ~ 6.25. (See the 
eigenvalues of the linearized Eliashberg equation for each 
channel A s and plotted in Fig. 6(d).) This suggests 
that the next-nearest-neighbor pairing of d x i_ y i orbital 
plays important role in stabilizing s-wave channel. This 
is a good example to show the convenience of our param- 
eterization in obtaining an intuitive picture for the gap 
function. 
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Fig. 6. (Color online) (a-c) Doping dependence of A33(fe). (d) 
Doping dependence of 03, 63 and C3. Eigenvalues of the linearized 
Eliashberg equation for s-wave pairing channel X s and for d-wave 
pairing channel are also plotted. 



In summary we have investigated the structure of the 
gap function in iron pnictides. We have shown that the 
combination of the orbital representation and real space 
view leads to a powerful and convenient expression for 
the gap functions. This expression is powerful enough to 
reproduce the results of microscopic calculation and to 
discuss the variety of the gap functions. At the same time 
it is simple enough to make an intuitive picture and to 
calculate physical quantities in the superconducting state 
through eq. (2). Calculation of the physical quantities 
in this formalism should give interesting information on 



whether the gap function is s-wave or d-wave, or the gap 

function is s ++ or s^ Similar analysis on 122- and 11- 

system will also give important information. 
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